Perturbative HFB model for many-body pairing correlations 
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We develop a perturbative model to treat the off-diagonal components in the Hartree-Fock- 
Bogoliubov (HFB) transformation matrix, which are neglected in the BCS approximation. Applying 
the perturbative model to a weakly bound nucleus 84 Ni, it is shown that the perturbative approach 
reproduces well the solutions of the HFB method both for the quasi-particle energies and the radial 
dependence of quasi-particle wave functions. We find that the non-resonant part of the continuum 
single-particle state can acquire an appreciable occupation probability when there exists a weakly 
bound state close to the Fermi surface. This result originates from the strong coupling between the 
continuum particle state and the weakly bound state, and is absent in the BCS approximation. The 
limitation of the BCS approximation is pointed out in comparison with the HFB and the present 
perturbative model. 

PACS numbers: 21.60.-n,21.60.Jz,21.10.Pc 



I. INTRODUCTION 

The pairing is the most important correlation beyond 
the independent particle approximation in the nuclear 
mean field P, Q . For the description of nuclei close to 
the beta stability line, the BCS approximation has been 
successfully employed in order to take into account the 
pair correlation. Although this approximation does not 
treat the interplay between the particle-hole and the pair- 
ing channels in a fully consistent manner , one can ex- 
pect that the defect is minor in well bound nuclei. As 
the Fermi energy approaches to zero, however, virtual 
scattering of nucleon pairs to the positive energy spectra 
may significantly alter the situation. If the continuum 
states acquire an appreciable occupation probability, in 
the BCS approximation, the particle density does not 
tend to zero asymptotically. In this case, the unphysi- 
cal particle gas surrounding the nucleus appears 0. 0.0. 
and one has to use the Hartree-Fock-Bogoliubov (HFB) 

method [lH0,Ellllllli, which is theoreticall y more 
robust than the BCS approximation. 

Recently, Sandulescu et al. have argued that the non- 
resonant part of the single particle continuum spectra can 
be neglected in the pairing correlation [l,[l(J[l3 (see also 
Refs. dill). Since the wave function for a resonance 
state is well localized inside the nucleus, the unphysical 
particle gas is negligible if the non-resonant part is ex- 
cluded from the positive energy space. By taking only a 
few narrow resonance states into account, the authors of 

Refs. HE [HE! have demonstrated that the BCS ap- 
proximation leads to similar results to the HFB method 
concerning the binding energy as well as the root mean 
square radius both for stable and weakly bound nuclei. 

Those results, however, do not necessarily mean that 
the BCS approximation provides a good representation 
for quasi-particle wave functions. In fact, Dobaczewski 
et al. have shown Q that the smaller part of the two 
component quasi-particle wave functions in the BCS ap- 



proximation behaves quite differently from that obtained 
with the HFB method, although the larger component is 
similar to each other. Since transitions are usually more 
sensitive to the details of wave function than energies, a 
more intriguing test of the BCS approximation should be 
for quasi-particle wave functions. 

The purpose of this paper is to discuss how the BCS 
approximation can be improved while retaining its intu- 
itive picture. To this end, we start from the full HFB 
equations and simplify them by the BCS approximation. 
We then propose a model which treats the effect beyond 
the BCS approximation by the perturbation theory. We 
perform detailed comparisons between the three models, 
the BCS, the HFB and the perturbative approximations 
for the quasi-particle wave function in weakly bound nu- 
clei, and clarify the limitation of the applicability of the 
BCS approximation. 

The paper is organized as follows. In Sec. II, we give 
the HFB equations in the coordinate space, and formu- 
late the perturbative approach to HFB. In Sec. Ill, we 
apply the perturbative approach to neutron states in the 
weakly bound nucleus 84 Ni, and discuss the applicability 
of the perturbation formulas. The summary of the paper 
is then given in Sec. IV. 



II. PERTURBATIVE APPROACH TO HFB 

Let us start with the HFJB equations in the coordinate 
space representation 
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where h and A(r) are the mean- field Hamiltonian and the 
pairing potential, respectively. Here, we have assumed 
that the nucleon-nucleon interaction is a zero range force 
so that the mean-field and the pairing potentials are lo- 
cal. A and E k are the Fermi energy and the quasi-particle 
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energy, respectively. In order to solve the HFB equations 
l[TJl. we expand the quasi-particle wave functions Uk(r) 
and Vfe(r) on the Hartree-Fock (HF) basis (fi, 



We next consider the effect of the off-diagonal compo- 
nents of Ay- . To this end, we rewrite the HFB Hamilto- 
nian Q as 
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where the wave function ipi satisfies htpi — eiipi. This 
leads to an eigenvalue problem of the matrix Hhfb which 
is given by 
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where 



Ay- = / drtf(r)A(r)<^(r). 



(-5) 



This is the so called two-basis method for the HFB 
equations E GS E3 

The BCS approximation is achieved by neglecting the 
off-diagonal components of Ay in Eq. 7]. The resul- 
tant equations, 



with 



(^BCs)ij 

(AW)y- 



Whfb = Wbcs + AH, 



(e, - A) An 
A„ ; (-a + A) 





(14) 



(15) 



Ay- - AuSij J ■ 111,1 



The eigen functions of Hbcs form the unperturbed basis, 
which is given by 
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Treating AH by the standard perturbation theory, we 
find that the correction to the quasi-particle wave func- 
tion and the quasi-particle energy is given by, 
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can be solved easily. The solutions are the well known 
BCS solutions [l|, 
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respectively, to the leading order. Here Mi is introduced 
to normalize the wave function \<pi). Notice that the 
first order correction to the quasi-particle energy van- 
ishes. The matrix element of ATI is given by, 



with 
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From Eqs. 0, ©, an d ©, it is apparent that the quasi- 
particle wave functions in the BCS approximation are 
given by 



U k (r) = «f cs ^(r), 
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Thus, the two components of the quasi-particle wave 
function are simply_proportional to each other in the BCS 
approximation 
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We call the perturbative model the p-HFB. 



III. APPLICABILITY OF THE 
PERTURBATION FORMULAS 

We now apply the perturbative approach, the p-HFB, 
to a drip-line nucleus in the 84 Ni region. Following Ref. 
0, we take a Woods-Saxon form for the mean-field and 
the pairing potentials, that is, 

V ws (r) = -V f(r) + V ls ~^-l-s (22) 

r ar 

(23) 
(24) 



A(r) = A /(r), 

m 



i 



1 + exp[(r - R )/a] 



We use the parameters, V =38.5 MeV, Vu=l4 MeV- fin, 
i?o= 5.63 fm, and a=0.66 fm in order to simulate the 
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FIG. 1: The matrix elements of the pairing potential A(r). 
These are plotted as a function of the single-particle energy 
€j for the state \<pj) with a fixed bound state (^»|, which 
has the same total and orbital angular momenta j and / as 
the state ifij. The open circles, the filled squares, and the 
open triangles are for the neutron lsi/j state at —31.6 MeV, 
the 3«i/2 state at —0.634 MeV, and the lgg/2 state at —4.11 
MeV, respectively. The diagonal components are denoted by 
the arrows. 



neutron potential around the 84 Ni region. The strength 
of the pairing potential is determined so that the average 
pairing gap defined as @ 



A 
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f™r 2 drf(r) 



(25) 



is equal to 1.0 MeV. The Fermi energy A is an input 
parameter in this simplified model, and we take A = —0.5 
MeV in the calculations shown below. For simplicity, the 
positive energy solutions of the mean-field Hamiltonian 
are treated by discretizing them with a box of 30 fm. 

We first discuss the matrix elements Ay of the pair- 
ing potential A(r) on the Hartree-Fock basis. Figure 1 
shows the matrix elements as a function of the single- 
particle energy of the state j for a fixed state i. For well 
bound states, such as the lsi/2 state at e = —31.6 MeV, 
they are all close to zero except for the diagonal compo- 
nent, which is denoted by the arrow (see the open circles). 
This justifies the success of the BCS approximation for 
stable nuclei. For weakly bound states, in contrast, the 
off-diagonal components are the same order as the diag- 
onal component. Specifically, the coupling between the 
weakly bound state and the positive energy states persist 
up to about 30 MeV. This is the case both for the s\/2 
state (the filled square) and for a state with larger angu- 
lar momentum such as the <?g/ 2 state (the open triangle). 
Since the BCS approximation neglects completely the off- 
diagonal components of the pairing potential, a careful 
investigation of the effects of the off-diagonal components 
is therefore necessary for weakly bound nuclei. 

Let us now solve the HFB equations by diagonalizing 
the HFB Hamiltonian J3J and compare with the BCS ap- 
proximation as well as with the p-HFB method discussed 



TABLE I: The quasi-particle energy E and the occupation 
probability v 2 for the four lowest neutron S1/2 quasi-particle 
states obtained with several methods. For the BCS method, 
the corresponding single particle level and the single particle 
energy are also tabulated. 





E (MeV) 


v J 


s.p. state 


e (MeV) 


HFB 


0.421 
1.045 
2.327 
4.259 


0.553 
0.039 
0.010 
0.003 






BCS 


0.462 


0.645 


3Si/ 2 


-0.634 




1.010 


0.007 


4si/ 2 


+0.509 




2.305 


0.004 


5S1/2 


+1.803 




4.244 


0.002 


6S1/2 


+3.742 


p-HFB 


0.400 
1.039 
2.318 
4.250 


0.606 
0.033 
0.005 
0.001 







in the previous section. We include the single-particle 
levels up to e— 50 MeV. We have checked that the re- 
sults do not change significantly even when higher lying 
states are included. Table 1 compares the HFB method 
with the BCS approximation for the quasi-particle en- 
ergy E and the occupation probability v 2 for the four 
lowest neutron Si/ 2 quasi-particle states. The occupa- 
tion probability is defined in terms of the quasi-particle 
wave function as 



dr\V k {r)\ 2 



(26) 



Also shown in the table is the corresponding single- 
particle state and its energy for the BCS approximation. 
We find that there is a good correspondence between the 
HFB and the BCS solutions, especially for the quasi- 
particle energy. One can thus assign unambiguously a 
particular single particle state to each HFB state as an 
unperturbative state. 

Figure 2 shows the difference in the quasi-particle en- 
ergy from the HFB energy for the neutron S1/2 states 
as a function of the energy of the corresponding single- 
particle state. The dashed and the solid lines are ob- 
tained with the BCS approximation and the perturba- 
tion theory, respectively. Surprisingly, the quasi-particle 
energy is off only by 0.04 MeV at most in the BCS ap- 
proximation, and the approximation seems to work well 
as long as the quasi-particle energy is concerned. This is 
in accordance with the finding in Refs. 0, 0, fill Il3| . 
The perturbative method (p-HFB) further improves the 
result, and the difference in the quasi-particle energy be- 
comes even smaller (see the solid line). Figure 3 com- 
pares the quasi-particle wave functions for the first three 
quasi-particle states. The left and the right panels show 
the V(r) and the U(r) functions, respectively. The solid, 
the dotted, and the dashed lines denote the results of the 
HFB method, the BCS approximation, and the p-HFB 
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FIG. 2: The deviation of the neutron S1/2 quasi-particle en- 
ergy, obtained with the BCS method (the dashed line) and 
the perturbative method p-HFB (the solid line), from the 
corresponding HFB energy as a function of the single-particle 
energy. 



FIG. 4: The occupation probability for the gg/2 states as a 
function of the corresponding HF energy e. The solid and 
the dotted lines are results of the HFB method and the BCS 
approximation, respectively, while the dashed line is obtained 
by the perturbative method, p-HFB. 
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FIG. 3: The quasi-particle wave functions for the neutron 
s-i/2 states obtained with the three methods, HFB, BCS, and 
p-HFB. The left and the right panels are for the lower (V) and 
the upper (U) components of the wave function, respectively. 
The solid line is the solution of the HFB method, while the 
dotted line denotes the wave function in the BCS approxima- 
tion. The result of the perturbative method p-HFB is shown 
by the dashed line. The top, middle, and bottom panels are 
for the 5si/2, 4si/2, and 3S]y 2 states, respectively. 



method, respectively. As noted by Dobaczewski et al. 3, 
the smaller component of the wave function in the BCS 
approximation is inconsistent with the HFB wave func- 
tion. Notice, however, that the perturbative method p- 
HFB dramatically improves the radial dependence of the 
wave function, suggesting that the BCS approximation 
still provides a good starting point to begin with. For 
the larger component, although the difference between 
the HFB and the BCS wave functions is small (but see 
below for the wave function for the lgg/2 state), the per- 
turbation again improves the shape of the wave function 
in a consistent way. The results of the p-HFB method 
are summarized in Table 1. 

We next discuss the quasi-particle states which have a 
relatively large angular momentum. As an example, Fig. 
4 compares the occupation probability for the gg / 2 states 
obtained by the HFB method (the solid line), the BCS 
approximation (the dotted line) and the perturbative ap- 
proach, p-HFB (the dashed line). Again, the p-HFB 
well reproduces the result of the HFB method. Notice 
that the occupation probability for the 3<7g « state at the 
HF energy of e =3.16 MeV is much larger than that in 
the BCS approximation, which essentially leads to a zero 
occupation probability for this state. The unperturbed 
quasi-particle energy, i.e., the quasi-particle energy in the 
BCS approximation, for this state (3g 9 / 2 ) is 3.66 MeV, 
which is close to the quasi-particle energy for the lgg/2 
state (3.73 MeV). Therefore, according to the perturba- 
tion theory, these two states strongly couple to each other 
by the off-diagonal components of the pairing potential. 
The wave functions for these two states are plotted in Fig. 
5. Indeed, the shape of the wave functions are similar to 
each other. In the HFB theory, hole states appear as 
narrow resonance states when the quasi-particle energy 

is larger than - AHHH0HI3. Width ° f ^ HFB 
resonance states of this kind originates from the coupling 
between the hole and the particle states due to the off- 
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FIG. 5: Same as Fig. 3, but for the l<?g/2 (the upper panel) 
and the 3<?9/2 (the lower panel) states. 

diagonal components of the pairing potential. Namely, 
the particle states which have a similar quasi-particle en- 
ergy as the hole state in the BCS approximation par- 
ticipate in the resonance state so that the quasi-particle 
wave functions behave similarly between the hole and the 
particle states within the resonance width|l8l. IliA l20| . 

This important mechanism of the particle-hole cou- 
plings is missing in the BCS approximation. Note that 
the 359/2 state is a non-resonant continuum state in the 
absence of the pairing. The resonance BCS method, 
therefore, is not expected to be adequate in the circum- 
stance where details of wave function plays an important 
role, for instance in the calculation of a strength func- 
tion with the quasi-particle random phase approximation 
(QRPA) IHI23. 

IV. CONCLUSIONS 

Within the mean-field approach, the pairing cor- 
relation is best incorporated with the Hartree-Fock- 
Bogoliubov (HFB) method. This method consistently 
takes into account the couplings between the particle- 



hole and the particle-particle channels, and thus theoret- 
ically robust. This method, however, is somewhat cum- 
bersome to solve, and the application of the fully contin- 
uum HFB calculation has been limited only to spherical 
systems so far. The HFB method, being based on the in- 
dependent quasi-particle approximation, also abandons 
a simple and intuitive single-particle picture. In this pa- 
per, we have developed a perturbative model to the HFB 
method, which provides an intuitive connection between 
the HFB and the BCS methods, and at the same time 
which can be solved much easier than the HFB. Apply- 
ing the perturbative approach to weakly bound nuclei, we 
have found that the lowest order perturbation reproduces 
well the results of the HFB method both for the quasi- 
particle energy and the radial dependence of wave func- 
tion. This suggests that the BCS approximation provides 
good unperturbative states, although it leads to incon- 
sistent quasi-particle wave functions with the HFB wave 
functions, especially for the smaller component. As good 
as the self-consistent HFB calculation, one may thus take 
into account the many-body pairing correlation by using 
the perturbative model on top of the BCS approximation 
after the convergence of the BCS calculation is achieved. 
This scheme might be easily applied to deformed nuclei 
for which the self-consitent HFB calculation is very duf- 
ficult. 

We have also pointed out that non-resonant scattering 
states in a mean-field potential can gain an apprecia- 
ble occupation probability when there is a weakly bound 
single particle state close to the Fermi energy. This im- 
portant coupling effect between particle and hole states 
cannot be incorporated in the BCS approximation. Ev- 
idently, the BCS approximation is inadequate for the 
QRPA calculations in weakly bound nuclei, where the 
details of the quasi-particle wave functions play an es- 
sential role. 
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